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(7^ We discuss neutrino fluxes and energy dependent flavor ratios of cosmic accelerators as a 

'""' function of the size of the acceleration region and the magnetic field, which are the pa- 

cn rameters of the Hillas plot. We assume that photohadronic interactions between Fermi 

, ^ accelerated protons and synchrotron photons from co-accelerated electrons (or positrons) 

O Isad to charged pion production. We include synchrotron cooling of the charged pions and 

^ successively produced muons, which then further decay into neutrinos, as well as the helic- 

ity dependence of the muon decays. Our photohadronic interaction model includes direct 

Q production, higher resonances, and high energy processes in order to model the neutrino 

O flavor and neutrino-antineutrino ratios with sufficient accuracy. Since we assume that the 

sources are optically thin to neutron interactions, we include the neutrino fluxes from neu- 

. !^ tron decays. We classify the Hillas plot into regions with different characteristic flavor ratios 

/\ over 20x24 orders of magnitude in R and B. In some examples with sizable magnetic flelds, 

^ we recover neutron beam, pion beam, and muon damped beam as a function of energy. 

However, we also flnd anomalous or new sources, such as muon beam sources with a flavor 

ratio v^w^Wj. of 1:1:0 in energy regions where the synchrotron-damped muons pile up. We 

also discuss the use of the Glashow resonance to identify p7 optically thin sources with a 

strong imbalance between vr"^ and tt^ production. We flnd that these can, in principle, be 

identifled in most cases in spite of the vr" contamination from high energy photohadronic 

processes and the mixing parameter uncertainties. 
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1 Introduction 



Neutrino telescopes |l]-[4j, such as IceCube or ANTARES, are designed to detect neutrinos 
from astrophysical sources. There are numerous candidate sources, the most prominent 
extragalactic ones being Active Galactic Nuclei (AGNs) (sl-fs] and Gammay Ray Bursts 



(GRBs) [9] , see Ref . 10 for a review and Ref . [TI] for the general theory of the astrophysical 
neutrino sources. For instance, in GRBs, photohadronic interactions are expected to lead to 
a significant flux of neutrinos in the fireball scenario. On the other hand, in AGN models, the 
neutral pions produced in these interactions may describe the second hump in the observed 
photon spectrum, depending on the dominance of synchrotron or inverse Compton cooling 
of the electrons. The protons in these models are typically assumed to be accelerated in the 
relativistic outflow together with electron and positrons by Fermi shock acceleration. The 
target photon field is often assumed to be the synchrotron photon field of the co-accelerated 
electrons and positrons. Furthermore, thermal photons from broad line regions or the 
accretion disc may serve as target photons, and there may be proton-proton interactions 
contributing to the neutrino flux. In this work, we focus on sources with a target photon 
field from synchrotron photons. Note that apart from AGNs and GRBs, there a many other 
candidates for neutrino sources. However, so far, no extraterrestrial high energy neutrino 
flux has been detected yet. This is, for sources optically thin to neutrons, consistent with 
generic bounds 



12 



13 which are just being touched by IceCube. Apart from that, there 
is also the possibihty to detect point sources for which the optical counterpart is absorbed, 
so-called "hidden sources" 



14 17 



We do not consider particular source properties, but we 
discuss generic point sources optically thin to neutrons. 

Apart from measuring muon tracks, neutrino telescopes may have some sensitivity to detect 



flavor 18,19 . The neutrinos are normally assumed to originate from pion decays, with a 
flavor ratio Uf.:v^:Vj. at the source of 1:2:0 (neutrinos and antineutrinos added) arising from 
the decays of both primary pions and secondary muons. However, it was pointed out in 
Refs. 11 20 that energy losses in strong magnetic fields, which dominantly affect the muons. 



change the flavor ratio at the source to 0:1:0. At low energies, neutron decays may dominate 
the flux, which leads to a flavor ratio 1:0:0. Therefore, one can expect a smooth transition 
from one type of source to the other as a function of the neutrino energy 21-23 , depending 



on the cooling processes of the intermediate muons, pions, and kaons. Recently, the use of 
flavor information has been especially proposed to extract some information on the particle 
physics properties of the neutrinos and the properties of the source; see Ref. 24 for a re- 
view. For instance, if the neutrino telescope has some flavor identification capability, such 
that it can distinguish between muon tracks and (electromagnetic or hadronic) showers, this 



property can be used to extract information on the decay ^,25 29] and oscillation 30 44 
parameters. Of course 



e.g. 



45,46 



the flavor ratios may be also used for source identification, see. 
Except from flavor identification, the differentiation between neutrinos and 
antineutrinos could be useful for the discrimination between p7 and pp induced neutrino 
fluxes, or for the test of neutrino properties (see, e.g., |27]). A useful observable may be 
the Glashow resonance process z/g + e~ — )■ W~ — )■ anything at around 6.3 PeV 18,34,47 



to distinguish between neutrinos and antineutrinos in the detector. Note that for photo- 
hadronic interactions, the vr"*" to 7r~ ratio determines the ratio between electron neutrinos 



and antineutrinos at the source. 

In the existing hterature, the flavor ratio is often taken to be a constant number independent 
of energy, with the above mentioned standard values taken as dogma. We demonstrate that 
even for "standard" sources, such as AGNs, flavor ratios in more realistic models could be 
anomalous in the sense that they do not necessarily fit the above picture. However, if the 
astrophysical parameters of the source, such as the extension of the acceleration region, 
are measured otherwise, such as by the variability timescale of the photon field flux, the 
energy dependence can be predicted in simple models. The idea of our model it to define 
the simplest possible toy model including a synchrotron target photon field and magnetic 
field effects on the secondaries. Naturally, it cannot describe all types of sources accurately. 
On the other hand, it does not rely on any astrophysical empirical relationships. Our main 
parameters are a universal injection index a for protons and electrons/positrons, the size of 
the acceleration region R and the magnetic field B. Our model is self-consistent in the sense 
that the target photon field is not put in by hand, but assumed to come from the synchrotron 
emission from electrons and positrons co-accelerated with the protons. Therefore, it relies 
on relatively few parameters (plus some assumptions). Pions are produced in photohadronic 
interactions between the protons and synchrotron photons, which leads to pion (and kaon) 
fluxes, and finally to neutrino fluxes. We include flavor and magnetic field effects, which 
means that we keep track of all flavors separately, include the magnetic field effects on pions, 
muons, and kaons, as well as the helicity dependence of the muon decays. The source is 
assumed to be optically thin to neutrons produced by photohadronic interactions, which 
means that we have to include the electron antineutrinos from neutron decays. Note that 
because we are only interested in the neutrino spectra, which are measured with too little 
statistics to exhibit time dependent features, we choose a time independent (steady state) 
approach. A similar self-consistent model can, for the specific case of microquasars, be 
found in Ref. ^8]. For other refined neutrino flux calculations, see, e.g.^ Refs. [7| [22|[49[[50] 

The parameters R and B can be directly related to the Hillas plot. In order to confine 
a particle in a magnetic field at the source, the Lamor radius has to be smaller than the 
extension of the acceleration region R. This can be translated into the Hillas condition for 



the maximal energy 51 



^max [GeV] - 0.03 -ri-Z-R [km] • fi [G] . (1) 

Here Z is the charge (number of unit charges) of the accelerated particle, B is the magnetic 
field in Gauss, and rj can be interpreted as an efficiency factor or linked to the characteristic 
velocity of the scattering centers. In this study, we consider the parameters R and B varied 
on the Hillas plot over 20x24 orders of magnitude. 

Since we study neutrino fluxes, flavor ratios, and neutrino-antineutrino ratios, we need a 
sufficiently fast model for the photomeson production which reproduces these features for 
power law-like spectra quite accurately. The often used A-resonance approximation 

._i_ j n + n^ 1/3 of all cases , . 

p + 7 ^ A ^ I ^ _^ ^0 2/3 of all cases ' ^^> 

is not sufficient for this purpose, since, for example, the vr" contamination is not included, 
and high energy processes affecting the spectral shapes are not treated at all. The best 



option may be the SOPHIA software 52 , a Monte Carlo simulation including not only the 
full final state kinematics, but also higher resonances, multi-pion production, and direct {t- 
channel) production of pions. This option is, however, too slow for our large scale parameter 
space scan, especially if smooth flavor ratio predictions are needed, and does not contain 
the helicity dependence of the muon decays, which is important for flavor ratio computa- 
tions (22). Therefore, we use Ref. [53] (model Sim-B) including the helicity dependent muon 



decays from Ref. 22 . This approach is largely based on the physics of SOPHIA. However, 
the kinematics are simplified, which can be seen mainly for sharp features in the proton 
and photon spectra at high energies, whereas for power law spectra, the treatment is very 
accurate. Compared to the approach in Ref. [54], which approximates the SOPHIA treat- 
ment analytically, the intermediate particles (pions, muons, kaons) are not integrated out, 
and magnetic field effects relevant for the flavor ratios can be included. For any given set 
of model parameters, our C-software NeuCosmA ( "Neutrinos from Cosmic Accelerators" ) 
can compute all secondary and neutrino flux spectra for one set of parameters in only a few 
seconds. This, of course, allows for applications which have not been possible in the past, 
such as the Hillas space scans in this study. In principle, our model can also be used as 
a (almost) online fit model for fluxes, because it interpolates among wide parameter space 
regions. 

This study is organized as follows: In Sec. |2} we describe our cosmic accelerator model, show 
a typical example, and discuss its limitations. The reader familiar with similar models may 



at least want to read Sec. 2J^, In Sec. [3| we discuss the source classification in the shock rest 



frame, i.e., without Lorentz boost, redshift, and flavor mixing. In Sec. |4| we then discuss 
the physics at the detector. We finally summarize in Sec. |5} 

2 The source model 

Here we discuss our model for the neutrino source. A (qualitative) model summary can 



be found in Sec. 2J^, whereas a more quantitative component description can be found in 
Sec. 2^ Then we discuss typical qualitative and quantitative results in Sec. 2^ In Sec. 2.4 
we comment on the limitations of our model. 

2.1 Short model summary 

We assume that protons and electrons (or positrons) are co-accelerated by Fermi processes in 
the source such that their injection spectra are power laws with the universal injection index 
a. The synchrotron radiation of the electrons in the magnetic field B of the source, assuming 
an isotropic pitch angle distribution, is the considered target photon spectrum for the pro- 
tons. Photohadronic interactions lead to charged pion, kaon, and neutron production. These 
particles decay further through weak decays into muons and neutrinos. Our photohadronic 
interaction model is based on Ref. [53] (model Sim-B) including direct production, different 
resonances, and high energy processes. The helicity dependent muon decays are taken into 



account as in Ref. 22 . The resulting neutrino spectra have to be Lorentz boosted (and 
redshifted) to obtain the spectra at the Earth. We assume that the source is optically thin 
to neutrons, which means that neutrons produced by the photohadronic interactions will 
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Figure 1: Flowchart describing the model. The functions Q{E) denote (injection) spectra per time frame 
[(GeVscm-^) ] and N{E) steady spectra [(GeVcm'^) ] derived from the balance between injection and 



losses or escape. Dashed arrows stand for solving the steady state differential equation Eq. ( 19 ) 



decay via beta decay and lead to an additional electron antineutrino flux. For all charged 
particles we include the effect of adiabatic and synchrotron energy losses, and we assume 
that these dominate such that energy losses due to inverse Compton processes, Bethe-Heitler 
pair production, and photohadronic interactions (the source is optically thin to neutrons) 
are sub-leading. Although not all spectral features of some sources may be described by 
this method, it has the advantage that the electron-positron (or proton) ratio is not needed 
as parameter, and that only the product of electron and proton injection spectrum enters 
the total normalization. 

To summarize, we show in Fig. [T] a flowchart describing the computation of the neutrino 
spectra in the shock rest frame (SRF). Here Q{E) (rectangular boxes) stand for (injection) 
spectra per time frame and N{E) (boxes with rounded corners) for steady spectra derived 
from the balance between injection and losses or escape. Dashed arrows stand for solving the 
steady state differential equation balancing energy losses and escape with injection. Weak 
decays and photohadronic interactions are denoted by the horizontal lines. The resulting 
neutrino spectra are marked by a double border and a colored background. Below the boxes 
we indicate in which decay the neutrinos are produced. 

The parameters of the source are all deflned in the SRF. We list them in Table [TJ together 
with units and typically used parameter values. The above mentioned magnetic fleld B 
is important for the synchrotron radiation and the related energy losses. The size of the 
acceleration region R plays an important role for the description of adiabatic losses and for 
the escape of particles leaving the source. The maximal energy of injected electrons and 
protons is limited where the energy loss rate exceeds the acceleration rate, which includes the 



Parameter 


Units 


Description 


Typical values used 


R 


km (kilometers) 


Size of acceleration region 


10ikm...l02ikm 


B 


G (Gauss) 


Magnetic field strength 


10-9 G... 10^5 G 


a 


1 


Universal injection index 


1.5. ..4 


-^e,mm 


GeV 


Minimal electron energy 


irie 


-^p,min 


GeV 


Minimal proton energy 


nip 


T] 


1 


Acceleration efficiency 


0.1 


r 


1 


Lorentz boost factor 


1 . . . 1000 


N 


a.u. 


Normalization/luminosity 


arbitrary units 



Table 1: Parameters of our model, units, description, and typical values used in this study. 



parameter r] (acceleration efficiency) from Eq. ([I]). The universal injection index a describes 
the index of the injection spectrum of protons and electrons. To fully describe these injection 
spectra the minimal energies of the proton (-Ep,min) and electron (-Ee.min) spectra and their 
normalizations have to be given. Our model is constructed such that only the product of 
these normalizations enters in the total normalization of the neutrino spectrum. Since we 
mainly discuss the spectra shapes and flavor ratios of the neutrino fluxes, we typically do 
not consider the normalization as separate parameter. The last parameter is a possible 
Lorentz factor F of a relativistic outflow which allows us to boost the neutrino spectra into 
the observer's frame at the Earth. We only consider this parameter in Sec. |4j The effects 
due to the cosmological redshift of the source are not explicitely discussed since they can 
be absorbed in a more generally defined Doppler factor (c/., Ref. [Tl]). 



2.2 Components of the model 

Here we quantitatively describe the key ingredients of our model. Unless noted otherwise, 
all quantities are given in the SRF. 

Synchrotron target photon field 



We calculate the synchrotron spectrum of the electrons in the Melrose-approximation 55 



averaged over the pitch angle. The power radiated per photon energy e by one particle with 
energy E, mass m and charge g in a magnetic field B is then given by: 



Ps{e,E) = l.t 



V^q^B 



1/3 



e ^/^" with Er 



3qB f E 



(3) 



16 Eomch \scj '' IQm \mc^ 

Considering a spectrum of particles radiating in a magnetic field we have to convolute with 
the spectrum of radiating particles as 



P{e) 



dEN{E)P,{e,E). 



The number of produced photons per time can be computed with: 



(4) 



(5) 



The steady photon spectrum, which is needed for the computation of photohadronic inter- 
actions, can be estimated by multiplying Q^ in units (GeVscm^) with the escape time 
tesc — R/c of the photons. The photon spectrum A^^ is then given by 

N,{e) = Q,{e) ^ (6) 

in units of (GeV cm^) . In the following, Q{E) will always denote spectra per time frame in 
units of (GeVscm^) , whereas N{E) will refer to steady spectra in units of (GeVcm^) 

Photohadronic interactions 

The accelerated protons interact with the synchrotron photons radiated by the electrons. 
In these interactions, pions, kaons, and neutrons are produced. We consider the following 
processes: 

P + 7 -)■ TT + p , (7) 

p + 7^i^+ + A/S. (8) 

Here p' is either a proton or neutron, tt are one or more neutral, positive, or negative charged 
pions, and A and S are resonances. Note that the effect of kaon decays is usually small. 
However, kaon decays may have interesting consequences for the neutrino flavor ratios at 
very high energies, in particular, if strong magnetic fields are present (the kaons are less 
sensitive to synchrotron losses because of their larger rest mass) l2ll|23]. Therefore, we 
consider the leading mode: K~^ production (for protons in the initial state) with the decay 
channel leading to highest energy neutrinosjj Note that at even higher energies, other 
processes, such as charmed meson production, may contribute as well. 

We follow the description of the photohadronic interactions in Ref. [53] (Sim-B), based 
on the physics of SOPHIA [52]. The treatment, taking into account resonant (including 
higher resonances), direct, and high energetic multi-pion production, allows for an efficient 
computation which is rather precise for power law-like spectra. It predicts well neutrino to 
antineutrino ratios and flavor ratios. The interaction between an arbitrary proton spectrum 
Np{Ep) and an arbitrary target photon field N^{e) is described by 

oo oo 

Q,(E,) = j ^ Np{Ep) j deN^{e)R,{x,y) (9) 

E,, ^th^P 

2Ep 

with X = Efj/Ep the fraction of energy going into the secondary, y = {Epe)/mp (directly 
related to the center of mass energy), the response function Rb{x,y), and eth the threshold 
for the photohadronic interactions. In our calculations, we use model Sim-B, which is 
characterized by the factorized response function 

R^ix, y) = 6{x - x'^) Ur F{y) with f^{y) = ^ / ^^^ ^r o^^^(e.) (10) 

eth 



"'^The contributions from K and K^ are about a factor of two suppressed 22 , and K'^ has a leading 
two body decay mode into neutrinos. 

6 



summed, for instance, over 23 different interaction types (IT) per pion species. For details. 



7r+ 


-> 


//+ -> e+ + z/g + z/^ , 


71^ 
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/i" ^ e" + z/e + z/^ , 


K+ 


-> 


/^^ + '^M , 
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-)■ 


p + e" + z/e . 



see Ref. 53 . With tlie factorization of the response function only a single integration is 
necessary to determine the spectra of the produced secondaries. This fact allows for a fast 
computation. From Eq. (|9|, it is also easy to see that in Qb{E) only the product of the 
normalizations of the proton and photon (radiated from electrons) fluxes enters. 

Weak decays 

We let the produced particle species, namely pions, kaons, and neutrons, decay weakly into 
neutrinos as 



:ir 



(12) 
(13) 
(14) 

Note again that for kaons, the branching ratio for this channel is about 64%. The second- 
most-important decay mode is K^ — )■ tt^ + tt^ (20.7%). The other decay modes account 
for 16%, no more than about 5% each. Since interesting effects can only be expected in 
the energy range with the most energetic neutrinos, we only use the direct decays from the 
leading mode. 

Since the decays of muons are helicity dependent, we distinguish between left and right 
handed muons (see Ref. [22]). In general, in case of ultra-relativistic parents of type a, the 
distribution of the daughter particle of type b takes a scaling form in order to obtain for the 
energy spectra 

summed over all parent species. The functions Fa^b for pion, kaon and helicity dependent 
muon decays can be read off from Ref. [22] (Sec. IV). For decays of neutrons produced by 
photohadronic interactions we use (c/., Ref. [22]): 

^n^Pe = ^ (I; - ^«) • (16) 

with x„ = 5.1 X lO^''. Since our source is assumed to be optically thin to neutrons and 
the neutron lifetime is relatively short, we assume that any produced neutron will decay 
between production and arrival at the Earth, either inside or outside the source. Thus, the 
neutrino spectrum can be directly related to the neutron injection spectrum as 

Q,^{E,) = —qJ^) . (17) 

Xn \ Xn J 
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Particle species 


Synchrotron Adiabatic 


Escape (region) Decay 


Electron/proton 
Target photons 
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X X 
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X 

X 



Table 2: Considered energy loss and escape timescales in the steady state for the considered particle 
species. 



Energy losses and escape 

For the charged species we consider possible energy losses, such as synchrotron radiation, 
before they decay. The kinetic equation for the particle spectrum, assuming continuous 
energy losses, is given by (see, e.g., Ref. J56J): 



with tcsc{E) the characteristic escape time, h{E) = —Et^J,^ with t^olsi^) ~ —^/EdE/dt the 
rate characterizing energy losses, Q{E) the particle injection rate [(GeVscm^) ] and N{E) 
the steady particle spectrum [(GeVcm^) ]. For different energy loss or decay mechanisms, 
the inverse timescales are summed over. Because the statistics in neutrino observations is 
not expected to be high enough for resolving time dependent features, we assume that energy 
losses and escape are in balance with the particle injection and that dN{E)/dt = 0, which 
leads to the so-called steady state spectrum N{E) described by the differential equation 

QiE) = ^{b{E)N{E)) + ^. (19) 

This steady state is another key ingredient for the efficient computation of neutrino fluxes, 
since time dependent calculations, which are the state-of-the-art for photon flux predictions, 
are extremely time intensive. Let us now discuss the special case of no energy losses, such 



as for photons, i.e., b{E) = 0. Then we can easily solve the differential equation in Eq. (19) 
in order to obtain 

iV(E)=g(E)tesc, (20) 

which is just the same equation as Eq. ^. 

We summarize the considered energy loss and decay timescales for the different particle 
species in Table 12} In the model, we consider the timescales for synchrotron (charged 
particles) and adiabatic (particles coupling to the plasma) energy losses and for escape 
from the region (neutral particles) and decays (unstable particles). The timescales can be 
approximated by the following formulas: 

^y'^^^'- 9 7reom4c5 ' ^'^ ^ R ' '' R ' ''^'^^ Etq ^ ' 

with g, E, m and tq the charge, energy, mass, and rest frame lifetime of the particle, B the 
magnetic fleld, and R the extension of the acceleration region. The expansion beta factor 

8 



in the comoving frame /5ex is approximated by one, which is a reasonable assumption for 
AGNs and GRBs (see Ref. pi). 



Injection spectrum and maximal energies 

The explicit form of our injection spectra for protons and electrons is: 

QpiEp) = 1^^ ^~" ^^''"'^ < ^P < ^P'-- (22) 

'^'=^''^^ = \0 else. ^ ^''^ 

Apart from the injection index a and the two model parameters -Ep,min and -Ee.min, the 
maximal energies Ep^^^x and -Ee.max have to be determined. We can estimate these energies 
by equalizing the characteristic timescale for acceleration of the particles with their energy 
loss timescale. For this computation, we use the acceleration timescale approximation fSl]^ 



rj, = v^^- (24) 

with r] an acceleration efficiency which depends on the acceleration mechanism, and Z the 
number of unit charges (see also Eq. ([T])). 

2.3 Typical results 

In this section, we first discuss what we can qualitatively learn about the secondary and 
neutrino spectra, and then we show a typical quantitative result. 

Qualitative shape of the neutrino spectra 

We have a spectral index a for the injection spectrum of the protons and electrons (in 
the following, the spectral index is always chosen positive). For large enough energies, the 
electron spectrum will be loss-steepened by one power due to synchrotron losses. Therefore, 
the index of the steady state spectrum of protons is a* = a and of electrons a^ = a + 1. 
Subsequently, we obtain an index a^ = (a* — l)/2 + 1 = a/2 + 1 for the photon injection 
spectrum Q^{E^). The spectral index of pions produced by photohadronics can be estimated 
by a^r = Op — a-y + 1. Hence, we obtain a^, = a/2. If the secondaries, such as pions and 
muons, do not loose energy, the daughter particles in weak decays have the same spectral 
index as the parent particles. Therefore, one obtains for the neutrino spectral index 

«. = f • (25) 



For a more refined discussion, see, e.g., Ref. 57 . Our assumption corresponds to the "efficient ac- 



celeration" regime therein. Although there could be additional constraints for the maximum energies for 
particles leaving the source from diffuse acceleration processes, for neutrino production, this discussion may 
only partially apply (the particles producing neutrinos interact before they leave the source). 

9 



o 



^e, max? ^— vJ-i 



15 



10 



10"- GeV 

10"' GeV 



NeuCosmA2010 



10' GeV 




O 
CQ 



^/?, max» ^— f .1 



NeuCosmA 2010 



10' GeV 



10' GeV 




Synchrotron cooling 
lO" GeV 



Figure 2: Maximal energy of electrons (left panel) and protons (right panel) in our model as a function 
of R and B for acceleration efficiency rj = 0.1. The solid (dashed) lines show the border between synchrotron 
loss and adiabatic loss determined maximal energy for 77 = 0.1(77 = 1). 



If there are synchrotron energy losses of secondaries, their (steady state) spectrum is loss- 
steepened by two powers above the critical energy 



Ec 



' 9 TT eo m^ c^ 



(26) 



which is defined as the energy where energy loss time and decay time are equal. This means 
that for energies higher than the critical energy the particles loose energy before decaying. 



Eq. (26) is valid for all charged secondaries with their corresponding masses, charges and 



lifetimes. Due to their shorter lifetime and their larger mass, the pions critical energy is 18 
times higher than for muons. The spectral index for neutrinos from pion a^ and muon a'^ 
decay respectively, is then given byjj 



a. 



K < e: 
2 E, > e: 



E, < E^ 

a:: = < f + 2 E^ > E^ 

f + 4 E, > E^ 



(27) 



The discussion of the spectral index for neutrino spectra applies to the energy range in 
which photohadronic interactions happen. From the maximal energy of the electrons follows 
the maximal energy of the photons ernax- This energy defines the minimal energy of the 
protons which still undergo photohadronic interactions, denoted by -E^ln' different to the 



•^For simplicity, we omit the fact that the neutrinos receive about one fourth (third) of the pion (muon) 
energy. 
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Figure 3: Inverse timescales for electrons and pions in the left panel and for pions, kaons and muons 
in the right panel. Here the parameters K = 10^'^ km, B — 10'^ G, a = 2, ii^p^min — "^p, £^e,min = "^e, and 
77 = 0.1 have been used. 



parameter -Ep^min describing the minimal energy of the proton spectrum (see Eq. (22)). Since 



the neutrinos obtain approximatively 1/20 of the proton energy, the index of the neutrino 
spectrum changes for energies below Ey < E^^jl^. The maximal neutrino energy can also 
be estimated by the maximal energy of the protons to -E^ax — -^p,max/20. 

In Fig. [2] we show the maximal energy of electrons (left panel) and protons (right panel) 
as a function of the radius R and magnetic field B for the acceleration efficiency r^ = 0.1. 
These energies are obtained from the balance between acceleration rate and loss rate. In 
the parameter region above the black solid lines synchrotron losses are dominant for the 
determination of the maximal energy, whereas below the lines adiabatic losses dominate. 
The dashed line corresponds to the solid line for r] = \. Changing in addition ri would 
change the maximal energy. In the region where synchrotron losses dominate the maximal 
energy, it is proportional to y^, and in the adiabatic regime it is proportional to 77. In 
general, it is clear that protons can be accelerated to higher energies than electrons due to 
their higher mass. Note again that it is well known that the maximum particle energies are 
limited by synchrotron losses in strong magnetic fields, see, e.g.^ Refs. 57,58 . However, 



the exact shape of the regions in Fig. |2] depends on the assumptions for the acceleration 
mechanism. 



A quantitative example 

In the following, we discuss the results for the model with the parameters B = 10^ G, 



R = lO^-^ km, a = 2, E, 



p,min 



m 



Ee,mm = "^e and T] = 0.1 (which is going to be test 
point 3 in the next section). In Fig. Isl we show the different inverse timescales (rates) for 
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electron and protons (left panel) and pions, kaons and muons (right panel) as a function 
of energy E. We can nicely see that the adiabatic cooling time is energy independent, the 
synchrotron loss rate proportional to E and the acceleration and decay rate proportional 
to 1/E. The dominant rate has always the largest values in Fig. [s] In the left panel we 
can read off the maximal energy for protons and electrons from the intersection point of the 
acceleration rate and the dominant energy loss rate (for comparison, see Fig.|2]). In the right 
panel, the critical energies for muons, pions and kaons can be read off from the intersection 
point of synchrotron loss and decay rate. Taking into account the maximal energy of the 
protons, we can infer from the critical energies that energy losses are for this parameter set 
not important for kaons, but relevant in the highest energy decade of the pion spectrum 
and the highest two energy decades of the muon spectrum. 

In Fig. Ill we show in the upper row the steady state spectra for pions (left) and muons 
(right) considering energy losses (solid curves) and without energy losses (dashed curves). 
Here the effect of energy losses sets in at around 10*^ GeV for pions and at 5 ■ 10^ GeV for 
muons, which are just the same values which can be read off from Fig. M (right panel). The 
muon energy spectrum changes the spectral index at both of these energies, because the 
parent pions are cooled above 10^ GeV. The qualitatively predicted behavior of the spectral 



indices leading to Eq. (27) can be clearly seen. In addition, there is a small "pile-up" effect 
in the muon spectra at around 10^ GeV, where the cooled muons from higher energies pile 
up to produce a flux higher than the uncooled one. 

In the lower row of Fig. |4l we illustrate the muon (left) and the electron (right) neutrino 
spectra coming from the different parents (c/.. Fig. [I]) comparing energy losses of the parents 
and no losses, where neutrinos and antineutrinos are added. For the muon neutrinos (left 
panel), one can clearly see that if there are large enough magnetic field effects, the different 
mass dependent cooling rates will lead to a spectral split of the corresponding neutrino 
spectra. The electron antineutrino spectrum from neutron decays is in this case sub-leading. 
It only contributes at low energies. In this specific example, the qualitative discussion above 
reproduces the spectral index of the neutrinos a,y = 1 below the cooling breaks, and the 
change of the spectra index at the breaks very well. However, the second break in the 
spectrum from muon decays is hardly visible. 

In the following we will, unless noted otherwise, fix some of the parameters [r] = 0.1, 
Ep^rain = f^p cLud -Ee,min = '"^e), siucc they do uot exhibit the most interesting effects. If 
one changes the acceleration efficiency, one changes the maximal energy of the protons and 
electrons and subsequently of all other particles from the decay chain. If one changes the 
minimal energy of the protons in the injection spectrum, the steady state spectrum has an 
index a^ = 1 below the minimal energy compared to a* = 2 above (if adiabatic energy 
losses dominate). This break translates into the pion and subsequently into the neutrino 
spectrum. In order to observe this effect in the important region of the neutrino spectrum, 
the Ep^j^in has to be chosen higher than about three orders of magnitude below the maximal 
proton energy, in our example, above 10^ GeV corresponding to a pre-acceleration 7^ ~ 10^, 
which is rather high. On the other hand, a larger value of -Ee.mm translates into a break in 
the steady state spectrum with a spectral index a^ = 2 below and a^ = 3 above. This break 
translates into a break in the photon spectrum around ec(-Se,min), see Eq. dsl), dependent on 
the magnetic field. The additional break will translate into the pion and neutrino spectrum 
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Figure 4: Pion (upper left), muon (upper right), niuon neutrino (lower left), and electron neutrinos 
(lower right) spectra for the model parameters B — 10'^ Q, R~ 10^-^ km, a = 2, i^p^min — "^p, £^e.inin = "^e, 
7/ ~ 0.1 (neutrinos and antineutrinos added). In the upper panels, the thick solid curves represent the total 
spectra considering energy losses, and the dashed curves neglect energy losses of pions and muons. In the 
lower panels, the thin solid curves show the results without energy losses (summed over all decay modes). 
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when Ep £c(-£'e,min) (with Ep < -Ep^max) is above the threshold for photohadronic interactions. 
The effect is then visible in the upper part of the pion and neutrino spectra. It turns out 
that the break in the neutrino spectrum has to be about three orders of magnitude below the 
maximal neutrino energy to be significant. For the model parameters chosen in this section, 
this translates into -Ee,mm ^ 1 GeV, corresponding to a pre-acceleration 7e ~ 2000. Of 
course, the pre-acceleration -Ee.min and -Ep^min depends on the specific model. For example, 
in some models a pre-boost comparable to the bulk Lorentz factor F may be a reasonable 
assumption. Since the effects on our results are, however, small and of secondary interest, 
we choose relatively low minimal injection energies. Note, however, that one could use a 
relatively high -Ee.min to model a (weak) break in a GRB photon field. 

2.4 Limitations of the model 

Our model is based on the simplest plausible set of assumptions, which means that it 
naturally cannot describe arbitrary astrophysical sources. For example, consider the target 
photon field. Since it is produced by the electrons with the same spectral index as the 
protons, the spectral index of the photon field is already fixed. Consider the often assumed 
photon field of GRBs with the spectral indices of 1 (below the break) and 2 (above the 
break), which is, in general, difficult to model with synchrotron radiation. In our model, 
a^ = 1 would require a hypothetical injection index for protons and electrons of a = 0, 



which again leads to a neutrino index of a^, = and not the typical a^ = 2 (see Sec. 2.3). 
To obtain the index a,y = 2 for the neutrino spectrum, we have to choose a = 4, but then 
the index of the photon spectrum is a^ = 3. We see that arbitrary combinations of indices 
of proton and photon spectrum are not possible since the injection index of electrons and 
protons is assumed to be equal and the photon index directly follows from the electron 
index. However, as mentioned above, one could use -Ee,min to model at least some spectral 
break. This problem arises mainly for GRBs and not for AGNs, since their photon fields 
can be described by synchrotron radiation much better. In that sense, our model is more a 
AGN-like model. However, to study the effects of secondary cooling, we use a = 4 in some 
cases, which may be unrealistically soft, but exhibits the features GRBs with respect to the 
pion spectra. 

In addition, with our model, we can not describe astrophysical environments where processes 
such as inverse Compton and Bethe Heitler pair production (with high radiation densities) 
are dominant or at least equivalent to the considered energy losses. However, in these 
cases, the assumption of optically thin sources to neutrons may also no longer be reasonable 
since the interaction rate of neutrons with photons depends on the photon field density. 
In high radiation densities, the interaction rate may be higher than the escape rate, which 
means that the neutrons would interact (and also cool via photohadronics) before leaving 
the source. For example, the synchrotron photon density is proportional to the absolute 
normalization of the electrons and the size of the source. Hence, in case of high radiation 
densities, one has to consider the normalizations of the injection spectra separately, or the 
ratio between protons and electrons, depending on the choice of parameters. 

In the calculation of the photon field we assume that we are in the slow cooling regime. This 
means, that the dynamical timescale of the system is shorter than the cooling timescale. 
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In the fast cooling regime, the dynamical timescale is longer than the radiative cooling 
timescale, as it may apply to the prompt emission of GRBs. The resulting photon spectrum 
is then softer by a factor of 1/2. The index of the photon spectrum is a!^^^ = al/2 + 1 



instead of aS, = {al — l)/2 + 1 (see 59 ). Note that we also neglect neutrino production 



through pp interactions. This process may contribute if high energetic protons run into cold 
protons, e.g., in GRBs during the afterglow. Including it into our model, however, would 
require more parameters. 

Finally, in our model we are exclusively interested in the target photon spectra for photo- 
hadronic interactions and not in the photon spectra observed at the Earth. To calculate 
the observed photon spectra, one has to take into account all radiation processes such as 
cascade radiation, synchrotron self absorption, and inverse Compton processes. However, 
note that this also implies that the observed photon field may not be the target photon field 
to be considered for photohadronic processes. Since our model is self-consistent, it does not 
rely on specific assumptions with respect to the multi-messenger connection. 

3 Neutrino fluxes and flavor ratios at the source 

Here we discuss the neutrino fluxes, flavor ratios, and (electron) neutrino-antineutrino ratios 
at the source. In particular, we discuss the classification of our sources as a function of R 
and B, as used in the Hillas plot. 

3.1 Source classification and the Hillas plot 

From Eq. ([I]), we can read off that if protons are accelerated {Z = 1) and the acceleration 
efficiency is high r^ ~ 1, the observation of the highest energetic cosmic rays with E ~ 
10^^ GeV leads to the necessary condition 

B[G]>l/3-10^%R[km])-^ (28) 

to produce such high energetic cosmic rays. Although we do not necessarily have to meet 
this condition for efficient neutrino production, it may serve as a hint for interesting source 
candidates. We show possible acceleration sites in the Hillas plot as a function of R and 



B in Fig. pi left panel (version adopted from M. Boratav; c/., Ref. 58 for a more detailed 
recent discussion). One complication in this type of figure is that R and B using the 
argument in Eq. ([I]) are potentially given in the SRF, whereas £'max is to be determined in 
the observer's frame. If the particles are accelerated in a relativistically moving environment, 
such as in a GRB fireball, this assignment is not trivial anymore, and the Lorentz boost 
r of the acceleration region must be taken into account by boosting R and B. We stick 
to the interpretation of R and B in the SRF, which means that (for rj = 1) the condition 



in Eq. (28) depends on the Lorentz boost of the source. This is illustrated by showing 
Eq. (28) for several selected Lorentz boosts in Fig. |5| left panel. We also include a GRB 
region "G" for illustration^] One can easily read off the figure that for the typical GRBs 
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GRB prompt emission; for typical parameters, which can be estimated from empirical relations, see, 



e.g., Refs. 60 61 . Note that the size of the acceleration region can be estimated from the timescale of the 



fluctuations in the GRB light curve (boosted into the SRF). 
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with r ^ 300, tlie maximuni energy in Eq. (28) can be reached, since the line T = 100 is 



partially exceeded. Note that some of the other regions also involve moderate boosts, which 
means that the interpretation in terms of R and B in the SRF has to be done carefully. 
For us, this interpretation is irrelevant, since we will use some points marked in Fig. [5] as 
test points with R and B given in the SRF, and we leave the interpretation to the reader. 



In addition, note that the condition in Eq. (28) is not necessarily matched by our model. 



Too large magnetic fields will lead to proton cooling by synchrotron radiation (such as for 
test points 1 and 2), which means that the maximum energy will be constrained by the 
synchrotron loss time in our model (c/.. Fig. [2]). In addition, we use r^ ~ 0.1 as acceleration 
efficiency in the following. 

The discussion at the source implies that all quantities are given in the SRF, such as 
magnetic field B and extension of the acceleration region R. In addition, the neutrino 
spectra are given at the source in arbitrary units (per energy, volume, and time) without 
Lorentz boost, redshift, or fiavor mixing. This discussion has the advantage that very few 
source specific parameters have to be taken into account. However, detector specifics, such 
as the detection threshold, cannot be taken into account for a source with arbitrary boost 
factor. For example, even if the neutrino energies are below the threshold at the source, the 
actual detectability depends on the boost factor, redshift, detection threshold, luminosity, 
and backgrounds. 

In the literature, typically the following three classes of neutrino sources are considered: 
Pion beams Neutrinos are produced by charged pion decays in the ratio VeWn'-Vr of 1:2:0 
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(neutrinos and antineutrinos added). This means that the electron to muon neutrino 
ratio Re = {Qu^+Que)/{Quf, +Qp^) — 1/2- In fact, corrections from the hehcity depen- 
dence of the muon decays lead to small deviations from this ratio, which, however, can 



be understood in terms of particle physics, as well as there can be pile-up effects 22 
The deviations depend on the input particle spectra (see, e.g., Ref. |53]) and are below 
the level of 10% in Rg. They are fully taken into account in our computations. 

Muon damped sources If the synchrotron loss time scale of the muons is shorter than 
the decay time scale, the muons loose energy before they decay. In this case, only 
neutrinos from pion decays are present, leading to u^'.i'^'.i'r of 0:1:0 (neutrinos and 
antineutrinos added), or R^ = {Qu, + Qu,)/{Qu^, + <5p,J — 0. 

Neutron beam sources Given the assumption of sources optically thin to neutrons, the 
neutrons leave the source and decay into protons, electrons, and electron antineutrinos. 
These neutrinos are typically found at low energies, unless the synchrotron losses of 
the pions and muons are so large {i.e., the magnetic fields are very large) that the 
neutrinos from neutron decays contribute significantly. In this case, we have Ve'-i^ii'-^r 
of 1:0:0 (neutrinos and antineutrinos added), or i?e = {Qu^ + Qu^)/{.Qu^ + Qv^,) -^ oo. 
Note that there are other potential sources of neutrons, such as neutrons produced by 



the photo-dissociation of heavy nuclei 62 , 63 . 

In addition, we find a new class of sources: 

Muon beams If the synchrotron loss time scale of the muons is shorter than the decay 
time scale, the muons loose energy before they decay. These muons may pile up at 
lower energies, where, similar to a neutrino factory, only neutrinos from muon decays 
are present. This leads to v^.v^w-r of 1:1:0 (neutrinos and antineutrinos added), or 
Re = (Qi/e + Qv^lKQv^ + Qp^) — 1- Typically, a muon beam comes together with 
a muon damped sources at higher energies: The muons missing at high energies are 
recovered at lower energies. Note that there may be other sources of muon beams, 
such as heavy flavor decay [D, Ds, B, etc.) dominated sources, where the pions and 
kaons interact before they can decay, whereas the heavy mesons have shorter lifetimes 
(see, e.g., Ref. j64| for a discussion). 

The above classification of sources in the literature is typically performed without consid- 
ering the energy dependence explicitely. In order to consider the full energy dependence, 
we have to take into account the energy ranges where most neutrinos will be detected, i.e., 
the energy ranges around the spectral peaks. The neutrino effective area roughly scales as 
E^Q (for muon tracks) below about 10 TeV even in the presence of Earth attenuation, as it 
can be seen in Fig. 6 of Ref. I65j. Therefore, we compute E^{Q^^ + QpJ and E'^{Q^^ + Qp^) 
separately, and then we determine the energy ranges where these fluxes are at least 1% of 
the respective maximum. The relevant energy range used for the classification of the source 
is then the union of these ranges. As we have explained above, we do not take into account 
the boost at this level. However, neutrino energies below 0.1 GeV at the source cannot 
contribute to the event rate, since these will be below the typical muon track thresholds 
(100 GeV for IceCube) even for extremely large F ~ 1000. Therefore, we do not consider 
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such low energies. Within the identified energy ranges, we then classify according to: If 
Vfr.v^w^ is between 0.9:0.1:0 and 1:0:0, i.e., R^ > 9, we have a neutron beam source. If 
Ue-.u^-Mr is between 0.1:0.9:0 and 0:1:0, i.e., Re < 1/9, we have a muon damped source 
source. If Re is between 0.45 and 0.55, we have a pion beam source. If Ve-v^Wr is between 
0.9:1.1:0 and 1.1:0.9:0, i.e., 0.9/1.1 < -Re < 1.1/0.9, we have a muon beam source. Other- 
wise, the source is mixed (undefined). We require that Re is within the respective range over 
at least one order of magnitude in energy, otherwise the source is again undefined. Note 
that this classification as a function of energy also admits sources with several well defined 
energy regions. For example, we find sources which are pion beams at lower energies and 
then translate into muon damped sources at higher energies. 

We show the result of this classification for a = 2 (injection index) in Fig. |5l right panel, as a 
function of R and B. Each of the points in this plot corresponds to a full model computation. 
The region in which no acceleration takes place or the photohadronic interactions are below 
threshold is marked by "no acceleration". Apart from that, all the characteristic sources 
mentioned above can be found. The classical pion beams cover most of the parameter 
regions, in particular, on the scales of galaxies and larger objects. Stronger magnetic fields 
lead to cooling, which means that the pion beams change into muon damped sources at 
higher energies (pion beam — )■ muon damped). The muon damped region consists of two 
parts: For small R < 10^ km, the muons loosing energies typically pile up at lower energies, 
leading even to separate "Muon beam — )■ muon damped" sources as a function of energy. 
For larger R, the increasing magnetic field at some point leads to the domination of muon 
damping around the peak independent of R. For even higher magnetic fields, both pions 
and muons strongly loose energy, which means that the neutrinos from neutron decays enter 
the relevant energy range, which eventually leads to neutron beams. 

In order to illustrate the classification of sources and typical examples even further, we 
show in Fig. [6] the muon neutrino (upper row) and electron neutrino (middle row) spectra, 
as well as electron to muon neutrino fiavor ratio (neutrinos and antineutrinos added, lower 
row) for several selected test points (in columns) from Fig. [s] In this figure, the energy 
ranges used for the source classification are gray-shaded according to the criteria above. 
Note that the strong jumps in the lower row above the energy classification windows may 
not be physical, since the fiuxes are vanishing there. In the upper two rows, we show the 
total muon or electron neutrino spectra, as well as the individual contributions from pion, 
muon, kaon, and neutron decays. In the lower row, the ratio leading to the classification 
as neutron beam, pion beam, muon beam, or muon damped source is also shown. Test 
point 11 is a typical example for a pion beam, as it can be seen in the lower middle panel. 
Here neutrinos are mostly produced by pion and muon decays, whereas kaon and neutron 
decays are mostly outside the classification range. The main spectral index of the neutrino 
spectrum is roughly one, as expected for a proton spectrum not dominated by synchrotron 
losses. The neutrino energies are moderately high (in fact, test point 11 has the lowest 
maximal energies of the shown pion beam test points). If the magnetic field is increased, 
such as in test point 13, muon damping leads to a clear hierarchy among the neutrinos 
from kaon, pion, and muon decays in terms of their energy; cf, upper right panel. This 
source is classified as muon damped above about 10^ GeV (lower right panel). Over the 
whole energy range, it qualifies as "pion beam — )■ muon damped" source, mainly because 
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the spectrum is relatively flat for a wide range in E^Q^. This flatness conies from the effect 
of proton synchrotron cooling. Compared to test point 13, test point 2 exhibits a much 
smaller acceleration region, which means that protons loose energy mostly by adiabatic 
cooling in our model. At lower energies, the muons damped from higher energies pile up, 
leading to the muon beam (c/., second panel in third row). The spectra split between 
neutrinos from muon and pion decays can be nicely seen in the second panel of the first 
row. Test point 12 is an interesting example of a mixed source (see fourth panel in third 
row). In this case, even the pions and kaons are strongly cooled (see fourth panel in first 
row), whereas the neutrinos from neutron decays also show up in the relevant energy range 
(see fourth panel in second row). The peak in the flavor ratio at about 10^ GeV corresponds 
to a narrow energy range where the neutron decays contribute, whereas for higher energies 
the muon neutrinos from kaon decays take over. Note that this effect from the kaon decays 
also leads to different conclusions than in Ref. [66] in the GRB-relevant range, where it 
shadows the neutron decays at high energies (including neutral kaons, the effect may be 
even stronger). In addition, note that in test point 12 the maximum neutrino energies are 
relatively small, and even smaller in test point 1, where synchrotron cooling does not allow 
for high proton energies. Although test point 1 may be below the detection threshold of 
a typical neutrino telescope, it is an interesting example of a source where the neutron 
decays dominate (lower left panel) because all charged particle species loose energy very 
quickly in the extremely strong magnetic field. Unlike other neutron beams discussed in the 
literature, where neutrinos are produced by the photo-dissociation of heavy nulcei, we have 
a genuine neutron beam here, which does not depend on the properties of the interacting 
particles other than being charged. However, it should be noted that neutral kaon decays 
may have an effect here (which we do not consider), and that the source will probably not 
be detectable anyway because of the too low energies. Test point 3, which lies close to three 
different regions (c/.. Fig. Isl), will be discussed in the next section separately. 



3.2 Dependence on injection index 

Let us now discuss the dependence of the source classification on the universal injection 
index a. While a = 2 may represent the often used standard case, many sources may have 
softer injection spectra, i.e., a > 2. Therefore, we show in Fig. [7] the source classification 
for four different values of a other than a = 2 in Fig. |5| a = 1.5, a = 2.5, a = 3, and a = 4. 
While the case a = 4 seems to be unrealistically soft at first, the pion and muon neutrino 
spectra in our model come then closest to the Waxman-Bahcall flux for GRBs [12] in the 
plateau region. Therefore, it is not surprising that the ORB region (c/.. Fig. [SJ left panel) 
lies within the "pion beam — )■ muon damped" region, as it is demonstrated in Ref. [20] . 

While some regions are almost unaffected by the injection index, such as the one around test 
points 5 to 11, there is some transition in other parts of the parameter space. Let us first of 
all discuss test point 3, which lies close to three different regions for a = 2 (c/.. Fig. [SJ), and 
for which the classification changes as a function of a. We show the neutrino spectra and 
flavor ratios for test point 3 and three different injection indices in Fig. |8} For a = 2, we have 
a muon damped source, because the lower energy range in the classification window exceeds 
the pion beam allowed range (see lower left panel). The reason is a combination of different 
factors, such as a small pile-up effect, the kinematics of the weak decays, and the non- 
20 
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Figure 7: Classification of sources in the Hillas plot for different injection indices as shown in the plot 
labels. 
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Figure 8: Muon neutrino (upper row) and electron neutrino (middle row) spectra, as well as electron 
to muon neutrino flavor ratio (neutrinos and antineutrinos added, lower row) for test point 3 for different 
a (in columns. The energy ranges used for the classification of the sources are gray-shaded. 
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negligible contribution from neutron decays (see middle left panel). For a = 3, the source 
qualifies as pion beam for most of the energy range. In this case, the contribution from 
neutron decays is relatively smaller. At high energies, however, we have a muon damped 
source (see lower middle panel). For a = 4, the muon damped range becomes shorter than 
the required order of magnitude (see lower right panel). Here the spectrum is flat in E'^Q^, 
which means that the statistics will be dominated by a wide energy range, and the muon 
damped part only gives a minor contribution. 

Apart from test point 3, there are a number of different interesting regions in Fig. [7] where 
the classification changes as a function of a. For example, consider test point 2, which is 
classified as "pion beam — )■ muon damped" for a = 4. For smaller a, however, the spectral 
shape changes, and, in the most extreme case a = 1.5, the neutrinos from muon and 
pion decays lead to two separated sharp peaks. In between these two extremes, the above 
mentioned combination of factors (as for test point 3) prevents the identification of the pion 
beam part, and leads to a muon damped source. On the upper end of the plot (large B), the 
dominance of neutron decays in the classification range depends on the injection index. For 
a > 2, the neutron beam region is pushed towards higher B until it vanishes, because then 
kaon decays contribute significantly. Another interesting region is the one around 10^^ km 
and 10^ G, which changes from "pion beam — )• muon damped" (small a < 2) over "pion 
beam" (medium a ~ 3) to an undefined source (large a > 4). Here the relative position of 
the peak of the spectrum (from the balance between injection and cooling) shifts to the left 
(to low energies) with increasing a, such that first of all the muon damped part leaves the 
classification energy range. For even larger a, strong pile-up effects and too low energies 
destroy the classification as pion beam. 

3.3 Optically thin p^ sources and the neutrino-antineutrino ratio 

Apart from the flavor ratios at the source, it is interesting to discuss which sources behave 
as optically thin (to neutrons) p'j sources. In order to identify such sources, one would 
use the Glashow resonance process at the detector, which is only sensitive to Ueij For an 
optically thin p7 source, protons interact with photons leading to mainly tt^, which then 



decay into /i^ and finally into Ue {cf., Eq. (11)). Therefore, almost no z/g are produced 



at the source. At the source, we therefore expect Qp^/Qu^ — for the "optically thin 
P7 source". Note that if the source was not optically thin to neutrons, ri'y (and other) 
interactions would break the asymmetry between Ue and Ug. Neutron decays would also 
increase Qp^/Qu^- In addition, if the source produced neutrinos by pp interactions, tt"*" and 
n^ would be produced in comparable amounts, which means that Qp^/Qu^ — 1. Therefore, 
we use the term "optically thin pj source" for Qp^/Qy^ — 0, because the source has to 
be both optically thin to neutrons, and most of the neutrinos have to be produced by p7 
interactions followed by charged pion decays. 

In practice, the above mentioned argument Qp^jQv^ — only holds in the A-resonance 
approximation Eq. ^ of the photohadronic interactions. If other processes, such as direct 



^Note that the Glashow resonance requires a particular neutrino energy, which, however, depends on the 
boost. Since we do not include a Lorentz boost yet, we classifiy the sources in the whole energy range in 
this section. 
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Figure 9: Classification of optically thin p'y sources (see main text) in the Hillas plot for different 
injection indices as shown in the plot labels. 



production and multi-pion production are included, the contribution from vr 
contamination from i^e] see, e.g., Refs. 52 
QuJQue ^ 0.2 as optically thin pj source 



leads to some 
Since we include these processes, we accept 



We use the same energy range for this classifi- 
cation as above for the flavor ratios, and we require that Qp^/Que ^ 0.2 over at least one 
order of magnitude in energy. Note that flavor mixing in combination with muon damping 
may introduce additional electron antineutrinos to the Qv^/Qu^ ratio, which we will discuss 
later. Without flavor mixing, however, any produced electron antineutrinos will be strongly 
suppressed. In this case, the classification as optically thin p7 source will be spoilt if neutron 
decays contribute significantly. 

We show the region of optically thin p7 sources for two different values of a in Fig. [9j Apart 
from large B, where neutron decays significantly contribute to the spectrum, and a region 
for a = 2 around test points 6, 8, 9, and 10, most of the sources qualify as optically thin p7 
sources. As we demonstrate below, for large B, the contribution of neutron decays prohibits 
the classification as such source. The undefined source region in the lower right corner of 
the left panel, however, comes from the significant contribution of high energy events in the 
photohadronic interactions, which break the 7r"^-7r~-asymmetry. There is only very little 
dependence on the injection index. However, this classification changes if flavor mixing is 
included within a slightly different analysis, as we will demonstrate later. 

We show several examples for selected test points and values of a in Fig. [T0| where the Ve and 
Uf. spectra are shown separately. For example, in the left column (test point 3 for a = 2) we 
find a typical optically thin p7 source within most of the classification energy range (gray- 
shaded). Here the neutrinos from neutron decays only dominant at the lower end of the 
classification range (see upper left panel). In the middle column, however, neutron decays 
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Figure 10: Electron neutrino and antineutrino fluxes (upper row), as well as electron antineutrino to 
neutrino ratio (lower row) for several selected test points (in columns) from Fig. ^ The energy ranges used 
for the classification of the sources are gray-shaded. 

contribute significantly over the whole energy range, which means that the source cannot 
be classified as optically thin p7 source. From the V(^ to z/e ratio, one can in this case not 
identify easily where the symmetry between electron neutrinos and antineutrinos actually 
comes from. The third column shows another example of a optically thin p7 source. In 
this case, however, neutron decays contribute significantly in a small energy window around 
10^ GeV, which leads to a "neutron spike" within the classification range. The regions where 
these neutron spikes occur are marked in Fig. [9j They can be much more prononounced 
than for the shown test point. 

4 Neutrino flavor ratios at the detector 

In this section, we consider the neutrino fiuxes at the surface of the Earth, including a 
possible boost of the acceleration region and neutrino fiavor mixing. We consider neutrino 
fiavor ratios as observables, which correspond to quantities in principle accessible to the 
detector. For these quantities, the absolute normalization cancels. We first show the impact 
of fiavor mixing on the observables. Then we illustrate how fiavor ratios can be used to 
extract information on the astrophysical parameters. And finally, we show the impact of 
Lorentz boost and fiavor mixing on the Glashow resonance. Note that we do not include 
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Earth attenuation effects in the present study. 

4.1 Impact of neutrino flavor mixing 

The easiest possibihty to measure fiavor ratios at the detector requires the identification 
of showers, which are much harder to detect than muons, and furthermore fiducial volume 
(not effective area) limited. Therefore, the statistics expected from showers is much smaller 
than from muon tracks. If we assume that electromagnetic (from Ue) and hadronic (from 
u-r) showers do not need to be distinguished, a useful observable is the ratio between muon 



tracks and showers 32 



R ^ 0Det + ^Det (29) 

where neutrinos and antineutrinos of each flavor are summed over. Here d^*^^ is the flux at 
the detector in units of cm~^ s~^ [sr~^], related to Q,^^ by a flavor independent normalization 
factor. Note, however, that neutral current events will also produce showers, which, in 
practice, have to be included as background. The benefit of this "flavor ratio" is that the 
luminosity of the source drops out. In addition, it represents the experimental measurement 
with the simplest possible assumptions. At IceCube, showers cannot be expected to be 
identified below about 1 TeV, which we consider as threshold for this observable. 

The neutrino flux at the detector is given by the neutrino flux per flavor (p^^^^^ at the surface 
of the Earth (we do not consider Earth attenuation effects here), modified by the effects of 
flavor mixing 

C = EEif^-n^/^^i''^'r'- (3°) 

In addition, we consider a possible Lorentz boost F of the acceleration region, which leads to 
a Lorentz boost of the neutrino energy E^ — )■ rii^j^jj Note that the transformation from the 
neutrino flux at the source to the flux at the detector includes a number of normalization 
factors coming from Lorentz boost, the integration over the production region, and the 
luminosity distance of the source. However, these will not affect the flavor ratios. 



Since in Eq. (30) the neutrino mixing angles appear, the uncertainties of the mixing angles 
will lead to uncertainties in the flavor ratios. In order to discuss these, we use the following 
values and current 3a ranges (see Ref. [67J): sin^ 6^23 = 0.5 {3a: 0.36 ... 0.67), sin^ 612 = 0.318 
{3a: 0.27 ... 0.38), sin^ 6'i3 = {3a: sin^ 6'i3 < 0.053). In addition, we consider future 
improved bounds (3o") on 613 from Daya Bay sin^ 613 < 0.012 and ^23 from T2K 0.426 < 
sin^^23 < 0.574, which represent the expected sensitivities at around 2015 |68|. Beyond 



these, a neutrino factory may improve the bounds even further: sin^ 6'i3 < 1.510 ^ 69 



and 0.46 < sin^ 6*23 < 0.54 70 {3a), which we label "2025". For 6*12, the errors may 



improve somewhat, but new generations of experiments are not assumed here. Since 6*12 is 



of secondary importance in Eq. (30), we do not consider such future improvements. 



^Often, a more general Doppler factor is used, which includes the Lorentz boost and the viewing angle 
and possibly redshift effects. The value of the Doppler factor D corresponds to the Lorentz boost factor F 
if the source is boosted at an angle l/F relative to the observer (redshift ignored). 
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Let us now discuss how the source identification changes in the presence of fiavor mixing. 
First of all, it does not make sense to consider an arbitrarily wide energy range, since 
the region close to the peak in E^Qy will contribute most significantly to the event rates. 
Therefore, we define, as in the previous section, the energy window which captures the 
upper two orders of magnitude in the fiux as "energy classification range". However, we 
use the muon neutrino fiux at the detector (neutrinos and antineutrinos added) for the 
determination of this window, since muon tracks are easiest to measure. Compared to the 
electron to muon neutrino ratio at the source (which may be very large for a neutron beam 
source), the observable R takes specific values for the considered source types and best-fit 
values: For the pion beam source, R = 0.5, for the muon damped source, R ~ 0.64, for the 
muon beam source, R ~ 0.44, and for the neutron beam source R ~ 0.28. We mark these 
values in the figures where applicable. 

We show in Fig. 
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the fiavor ratio R after fiavor mixing as a function of energy for test 
points 2, 11, and 13 from Fig. [s] (right panel). The solid curves show the result for the 
current best-fit values, whereas the shaded regions show the ranges if the neutrino mixing 
parameters are varied within their current or future (for 2015 and 2025) 3a ranges. Shaded 
rectangles mark the classification energy ranges where most events are expected. In addition, 
the assumed IceCube threshold is shown. The fiavor ratios for pion beams, muon damped 
sources, and neutron beams are marked. Here a Lorentz boost F = 1 is assumed, larger 
Lorentz factors would shift the spectrum (including classification range) to the right, but 
the IceCube threshold would remain fixed. 



Comparing Fig. 11 to Fig. [6] for the current best-fit values (c/., solid curves in Fig. 11), 
the energy dependence at the source translates into the energy dependence with the same 
features at the detector (only the vertical direction is inverted). Therefore, in principle, 
the source classes can be identified as a function of energy. The upper panel corresponds 
to a muon damped source, the middle panel to a pion beam source, and the lower panel 
to a pion beam — )■ muon damped source. If, however, the current parameter uncertainties 
are included in the considerations, a particular value of R cannot be assigned to a specific 
source type anymore. For example, consider the lower panel: If i? ~ 0.5, the source may be 
a pion beam, a muon damped source, or anything in between. As it can be read off from the 
figure, already the improved precision on the mixing parameters in about 2015 allows for a 
classification into the different source types. Even the muon beam range in the upper panel, 
R = 0.5 (pion beam) could be, in principle, excluded. Note that we have chosen 3a ranges 
in this figure, which is already quite conservative. In addition, we have chosen ^13 = as 
the best-fit, because otherwise we would need to discuss the measurement of ^cp as well. If 
Nature has chosen ^13 > 0, a neutrino factory may eventually be necessary to constrain 613 
and 6cp simultaneously. The main uncertainty, however, comes from 623- A different value 
of r would not change our results qualitatively, but instead shift a larger region beyond the 
assumed IceCube threshold. 

4.2 Measurement of astrophysical quantities 

Once neutrino fiavor ratios are measured in neutrino telescopes, it is conceivable to use such 
observations to constrain the parameters of the source. We know, for instance, that the 
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Figure 11: Flavor ratio R after flavor mixing as a function of energy for test points 2, 11, and 13 from 
Fig. [5] (right panel). The solid curves show the result for the current best-fit values, whereas the shaded 
regions show the ranges if the neutrino mixing parameters are varied within their current or future (two 
different assumptions) 3cr ranges (see main text for values) . Shaded rectangles mark the classification energy 
ranges where most events are expected (see main text). The flavor ratios for pion beams, muon damped 
sources, and neutron beams are marked. Lorentz boost factor F = 1 assumed. 
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Figure 12: Contour regions of R in the plane of energy versus magnetic field for a = 2 and R = 10^'^ 
km. The best fit values of the neutrino mixing angles were used. The region between the two dashed white 
curves is the energy classification range (see text). In the black regions (right corners) the neutrino fluxes 
go to zero so R is not well-defined. 

observation of supernova neutrinos is often used to constrain the physical processes that 
take place within the supernova, so the same idea can be applied also to the neutrinos 
produced in cosmic accelerators. 

Assume, for instance, that an optical (or gamma ray) counterpart observation of the source 
has allowed the measurement of F, R and a, e.g., via the time dependence of the photon 
spectrum. Can B be extracted from the energy dependence of the flavor ratio Rl In general, 
it is difficult to do so because R is not very sensitive to the magnetic field. In fact, we see 
from Fig. |5J that, at a given i?, the magnetic field can generally vary over several orders 
of magnitude without changing the classification of the region. But, close to a boundary 
between two regions, the magnetic field determines the classification of the source and, 
therefore, the expected value of R. It is close to such boundaries that the measurement of 
R may allow to constrain the magnetic field. 

Fig. 12 illustrates the variation of R with the energy and the magnetic field for a = 2, 
R = 10^'^ km (corresponding to TP3), and the best fit values of the neutrino mixing angles. 
Four regions are explicitly shown: R < 0.35, 0.35 < R < 0.45, 0.45 < R < 0.55 and 
R > 0.55. In addition, the energy classification range is shown by the dashed curves. From 
the figure we see that for energies around 1 TeV, R is typically larger than 0.45 independently 
oi B. It is also observed that the region i? > 0.55 is characterized by large magnetic fields 
in the energy range of interest for neutrino telescopes. Thus, if R were measured to be 
smaller than 0.55 for energies between 1 and 10 TeV, we could deduce, from this figure, 
that B < 10^ G. Clearly, it is possible to obtain some information about astrophysical 
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Figure 13: R a.s a function of the energy for three different values of the magnetic field. The bands 
indicate the expected uncertainty in 2015 (left) and in 2025 (right). The vertical lines mark the energy 
classification windows (see main text). 

parameters from flavor ratio measurements. 

It is important to analyze if this conclusion still holds when the uncertainties in the neutrino 
mixing angles are taken into account. To that end, we show in Fig 
the energy for three different values of B. As before we take a 
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i? as a function of 
2 and i? = 10^-^ km, but 



we now include the expected uncertainty in neutrino mixing parameters for 2015 (left) and 
2025 (right). From the figures, we notice that by measuring R over an energy range above 
1 TeV, one could indeed gain some information about the magnetic field, li B = 10^ G 
(upper shaded/red region), for instance, one expects a transition from a pion beam source 
to a muon damped source for energies between 10 and 100 TeV. But if S = 10 G (lower 
shaded/blue region), no transition is expected and R is never larger than about 0.5. We 
see then that the expected uncertainties in the mixing angles do not spoil the possibility of 
using R to obtain information on the magnetic field of cosmic accelerators. However, the 
information from 2015 is on the borderline on being useful, since the regions, which come 
from very different magnetic fields, somewhat overlap. 

4.3 Optically thin pj sources? 

In order to test the Glashow resonance at about 6.3 PeV, which is sensitive to electron 
antineutrinos only, it is useful to consider the flavor ratio 



LDet 



T 



ADet 



+ 



LDet 



(31) 



at the detector, representing the number of Glashow resonance versus muon track events. 
The fluxes at the detector, of course, include flavor mixing, which depends on the mixing 
parameters. Therefore, we need to include the uncertainties on the mixing parameters. 

The purpose of this process can be the discrimination between p7 optically thin sources, 
which produce more tt"*" than 7r~, and pp or optically thick sources, which produce compa- 
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Pion beam Muon damped Muon beam 



T for best-fit values 

Ideal pp source 0.50 0.28 0.64 

Ideal p7 source 0.22 0.00 0.36 

P7 source with 20% vr" contamination 0.31 0.09 0.51 
2^mm ^Qj, ideal pp source (marginalized over oscillation parameter ranges) 

Current bounds 0.38 0.15 0.56 

2025 0.47 0.23 0.62 



Table 3: Values of flavor ratio T used for "p7 opticafly thin" analysis (see main text for details). 

rable amounts of 7r+ and vr^ . Even after flavor mixing, the value of T is much smaller for a 
P7 optically thin source than any of the other mentioned sources. For instance, for the ideal 
muon damped p7 source, T ~ because no antineutrinos (of either flavor) are produced at 
the source. In order to reject the hypothesis of a source symmetric in vr^-vr^ production, we 
compute the theoretical minimal possible value T™'"^ for the 7r+-7r~ symmetric source after 
flavor mixing for a pion beam, muon damped source, and muon beam source separately, 
assuming that tt"*" and vr" are produced in equal amounts at the source. For that we vary 
the oscillation parameters in their 3a allowed ranges, or use the bounds from 2025 (these 
values correspond to the lower end of the marginalized T). The resulting values can be read 
off from Table [3| In the table, we also show the T for the used best-flt values for an ideal 
pp source (vr^-vr" symmetric) and ideal py source (Tr"*" only) expected from the A-resonance 
approximation only (c/., Eq. (^), as well as for the a pj source with a 20% tc~ contamina- 
tion, as it can be expected from high energy processes contributing to the photohadronic 
interactions. 

In order to identify an optically thin p7 source, we require T < T™™ at the detector for 
E = 6.3 PeV. Note that we do not include statistics and systematics in this discussion, we 
only discuss the theoretical possibility to use the Glashow resonance. From Table |3| we can 
read off that even for the current bounds within each source type, T < T™™ even for the 
py sources with a 20% contamination of vr" (compare the values from the third and fourth 
rows fllled with numbers). Therefore, a source identiflcation is, in principle, theoretically 
possible. This possibility is so far based on the assumption that the source can be identifled 
and that the flavor ratio at the source is "ideal" for each speciflc source. It is obvious from 
the values for T™™ that the identiflcation of the source is mandatory at least for pion beam 
and muon beam sources. For example, assume that T c^ 0.3 is measured without source 
identiflcation. In this case, it can be interpreted as muon damped pp, or pion beam pj 
source. Therefore, we require that the type of source (pion beam, muon damped or muon 
beam) can be identifled from the flavor ratio with the same energy ranges (at the source) 
as discussed in Sec. |3| over at least one order of magnitude. However, we choose the muon 
neutrino flux at the detector to determine the relevant energy classiflcation range. 

For each source type within its allowed energy range in ii^source, we compute the Lorentz boost 
factor range from P ^ 6.3PeV/-Esource for which the Glashow resonance is matched and for 
which T < T™™, i.e., the optically thin p'j source can be identifled. Note that compared to 
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Sec. 3.3, we focus on one particular energy (at the detector) for the p7 classification. This 
particular energy, however, corresponds to a range F™" ^ T ^ pmax ^^^ which a particular 
source can be used. Finally, we require 1 < F < 1000 for reasonable sources. Note that 
neutron beams are not useful for this type of analysis, since they are not produced by pj 
interactions. However, whatever effect destroys the 7i~^-7i~ asymmetry, makes T increase, 
i.e., acts in the same direction. 



We show the result of this analysis for a = 2 in Fig. 14 The left column shows the minimal 
allowable F and the right plot the maximal allowable F within the analysis range. The 
upper row is for the current bounds, the lower row for 2025. If there are several possible 
regions which can be used for the Glashow resonance (such as in the pion beam — > muon 
damped case), we choose the smallest and highest F from these regions. This figure (upper 
panels) is to be interpreted as follows: For example, for test point 6, the Glashow resonance 
can be used for 1 ^ F < 100, which is a reasonable range since these sources will be found 
around F ~ 1. The same argument applies to test points 4 to 10 for different maximal 
allowed values of F, which are the typical pion beam sources which can be used for this type 
of analysis, and test point 2 as muon beam — )■ muon damped source. Note that for test 
point 4, for which the Lorentz boost is often larger than 10, values as large as 1000 are also 
fine. Since the neutrino energies for test points 11 are lower, higher boost factors F > 10 
would be required, which might be unrealistic for this type of source. 



In Fig. 14, there are several gray-shaded regions marked "no useful source" where the 
Glashow resonance cannot be used. The reason is different in different parts of the parameter 
space. For large values of S, we have neutron beam sources or mixed sources, for which the 
discrimination between pp and p7 interactions is not applicable. Around B ~ 10^ G, the 
synchrotron cooling prevents neutrino energies high enough, which means that F ^ 1000 
would be required. The same applies to the region around R ~ 10^ km and B ~ 10^ G. 
Around test point 3 and at the spike at i? ~ 10^^ km and B ~ 10^'^ G the energies in the 
identified ranges are typically too high to match the Glashow resonance. Only for improved 
bounds on the mixing parameters (lower panels), a muon beam range at lower energies 
can be found in some part of the parameter space, where the source can be identified and 
the Glashow resonance is matched together with T < T™™ in a certain F range. For test 
point 3, the required F factors are still very high. Note that the energy classification window 



is slightly different than in Sec. 3.3 here. 



5 Summary and conclusions 

We have considered neutrino production by photohadronic interactions in cosmic accelera- 
tors. Our model assumes that protons interact with photons from co-accelerated electrons 
and positrons. We have included magnetic field effects, such as synchrotron cooling of the 
charged particle species, and flavor effects, such as the helicity dependent muon decays af- 
fecting the flavor ratios at the source and flavor mixing. In addition, we have considered 
neutrinos from kaons and neutrons produced in the photohadronic interactions. Because the 
spectral shapes, flavor ratios, and neutrino-antineutrino ratios, which we have investigated, 
strongly depend on the initial charged pion distributions coming from the photohadronic 
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Figure 14: Classification of optically thin p7 sources from the flavor ratio T in the Hillas plot for a — 2. 
The upper row is for the current bounds of the neutrino mixing parameters, the lower row for the 2025 case. 
The left column shows the minimal F of the source allowed for matching the Glashow resonance energy, the 
right column the maximal F (see main text for details). 
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interactions, we have used an efficient method to compute the secondary particle spectra 



based on Ref. 53 parameterizing the physics of SOPHIA. This approach includes direct 
production and high energy processes and allows for a computation of one model within a 
few seconds. We have used this approach for a large scale parameter space scan as a function 
of R (size of the acceleration region) and B (magnetic field) to study different sources as a 
function of the Hillas plot parameters. Another important parameter of our model is the 
(universal) injection index a. 

One of the quantities of interest has been the ratio between electron and muon neutrinos 
("flavor ratio") at the source, which characterizes its flavor composition before propagation 
to the Earth. For example, "pion beam sources" produce neutrinos in the flavor ratio 
Ue'Mf^-.Ur of 1:2:0, "muon damped sources" in the ratio 0:1:0 (the muons loose energy in 
the magnetic field before they decay), and "neutron beam sources" in the ratio 1:0:0 (the 
neutrinos are produced by beta decays). As a new type of source, we have found "muon 
beam sources" with the ratio 1:1:0, similar to a neutrino factory. In this case, the muons 
from higher energies pile up in the lower energy part of the spectrum, where pion decays 
hardly contribute. We have illustrated that a discussion of the flavor ratio only makes sense 
if one considers the corresponding spectrum simultaneously, since the region around the 
spectral peak (in E'^Q^y) will contribute most to the event rates. Therefore, we have focused 
on an energy range around the peak for the source classification, and we have required a 
unique classification over at least one order of magnitude in energy. For a = 2, our main 
result can be found in Fig. [5] Interpreting the test points in this figure in terms of the 
sources listed there, most sources on galactic scales perform as the classical pion beam. 
Neutron stars, on the other hand, are typical neutron beams because the strong magnetic 
field damp all charged particle species. However, the neutrino energies are probably too 
low to be detected. The most interesting regions may be around white dwarfs (TP 2) and 
active galactic nuclei (TP 3). For example, active galactic nuclei (TP 3) are classified as 
muon damped sources for a = 2 because neutron decays and magnetic field effects prohibit 
the classification of the pion beam energy range. For a = 3, however, they show both 
pion beam and muon damped ranges, and for even softer sources they perform as pion 
beams. The white dwarfs (TP 2), however, show for a = 2 a muon beam — )■ muon damped 
behavior, i.e., the muon neutrino spectra exhibits a spectral split coming from the muon 
damping. Although the target photon field of GRBs cannot be described within our model, 
we have shown that for pion spectra characteristic in GRBs the characteristic pion beam — )■ 



muon damped energy dependence from Ref. 20 is reproduced. Apart from the discussed 



flavor ratios, fully anomalous flavor ratios are found in some regions of the parameter space 
especially in regions where strong magnetic field effects compete with neutron decays (the 
neutrons are not affected by the magnetic fields). 

While the above classification has been performed at the source, the situation at the detec- 
tor is different. First of all, flavor mixing enters as a propagation effect. And second, the 
acceleration region may be Lorentz boosted with respect to the observer blue-shifting the 
neutrino energies, whereas the cosmological expansion redshifts the energies. The absolute 
energies at the detector are, however, relevant, since they have to be above the detection 
threshold. We have demonstrated that these effects do not, in principle, prohibit the de- 
tection of flavor ratios including the source classification for many characteristic test cases. 
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However, the knowledge on the neutrino mixing parameters can only be expected to be good 
enough after the T2K and Daya Bay results for serious applications. On the other hand, 
we have shown for a luminous point source, that some information on the source may be 
extracted from the flavor ratios. For example, if all the other parameters can be estimated 
from related measurements, such as the size of the acceleration region by the variability 
timescale in gamma ray observations, the flavor ratio allows in principle for independent 
constraints on the magnetic fields of the source. 

Finally, we have discussed if the Glashow resonance at 6.3 PeV, which is sensitive to elec- 
tron antineutrinos only and therefore to the neutrino-antineutrino ratio, can be used to 
learn something about the neutrino production mechanism. For example, in jry interactions 
significantly more vr"'" than 7r~ are expected to be produced in sources optically thin to 
neutrons. This asymmetry can be destroyed by subsequent neutron interactions in optically 
thick sources, a significant contribution of the high energy photohadronic processes, a signif- 
icant contribution of pp interactions, or neutron decays leading to a z/g typically not present 
in such sources. We have performed an analysis at both the source and at the detector. 
While at the source, our "p7 optically thin" condition is satisfied for most accelerators at 
least in some energy range, at the detector the Glashow energy has to be exactly matched, 
and the type of sources has to be identified to disentangle the effects from flavor mixing. 
This leads to constraints to the possible boost factors F of the sources. For example, the 
active galactic nuclei (TP 3) typically cannot be identifled as p7 optically thin. Sources on 
galactic scales or larger typically can be easily classifled. 

We conclude that flavor and neutrino-antineutrino ratios of astrophysical sources may be 
interesting to study both particle physics and astrophysics properties of the sources. Flavor 
ratios are an important possibility to test the magnetic flelds in the source, because these 
lead to a spectral split of the neutrino spectra from pion, muon, and kaon decays. On the 
other hand the neutrino-antineutrino ratio can be used to study the neutrino production 
mechanism. Propagation effects between source and detector do, in principle, not void these 
effects. On the other hand, they may be used to study non-standard neutrino properties, 
such as decoherence or decays. Therefore, the study of flavor ratios and the Glashow 
resonance at neutrino telescopes should be higher valued, and future instruments should be 
optimized for that once cosmic neutrinos are detected. 
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